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Magnetic fields correlated on several kiloparsec scales are seen in spiral galaxies. Their origin 
could be due to amplification of a small seed field by a turbulent galactic dynamo. We critically 
review the current status of the galactic dynamo, especially some of its problems and possible 
solutions. We also comment on the nature of seed magnetic fields, needed to prime the dynamo. 
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, I. INTRODUCTION 
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Magnetic fields in spiral galaxies have strengths of order few 10 _6 G, and are coherent on scales of several kpc [Q. 
In several disk galaxies, like M51 and NGC 6946, they are also highly correlated (or anti-correlated) with the optical 
»_i i spiral arms. How do such ordered, large-scale fields arise? One possibility is dynamo amplification of a weak but 
' nonzero seed field ~ 10~ 19 — 10 _23 G, provided the galactic dynamo can operate efficiently to exponentiate the field 
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by a factor ~ 30 — 40. We critically review here some of the issues relevant to the operation of the galactic dynamo, 
' the problems which arise and possible solutions. We also touch upon the origin of the seed magnetic field, needed for 
dynamo amplification. 

t-H ■ The evolution of the magnetic field, in the MHD approximation, is described by the induction equation 

>' 

Cn| i — = Vx (vxB- ? yVxB), (1) 

provided one assumes the usual form of Ohms law and neglects the displacement current term in Maxwells equation. 
Here B is the magnetic field, v the velocity of the fluid and r\ the resistivity. If ?y — > the magnetic flux through any 
area in the fluid is conserved during the motion of the fluid. The presence of a finite resistivity allows for a violation 
of such 'flux freezing' and the magnetic Reynolds number (MRN) R m — vL/tj measures the relative importance of 
' flux freezing versus resistive diffusion. Here v and L are typical velocity and length scales of the fluid motions. In 
most astrophysical contexts flux freezing greatly dominates over diffusion with R rn >> 1. 

B — is a perfectly valid solution of the induction equation. So there would be no magnetic field generated if one 
were to start with a zero magnetic field. The universe probably did not start with an initial magnetic field. One 
therefore needs some way of violating the induction equation and produce a cosmic battery effect, to drive curents 
from a state with initially no current. There are a number of such battery mechanisms which have been suggested 
[gj. All of them lead to only small fields, much smaller than the galactic fields. Therefore dynamo action, due to a 
velocity field acting to exponentiate small seed fields efficiently, is needed to explain observed fields. We first briefly 
comment on a cosmic battery before discussing dynamos in detail. 



o> 



C3 



II. COSMIC BATTERIES AND SEED FIELDS FOR THE DYNAMO 



The basic problem any battery has to address is how to produce finite currents from zero currents? Most astro- 
physical mechanisms use the fact that positively and negatively charged particles in a charge-neutral universe, do 
not have identical properties. For example if one considered a gas of ionised hydrogen, then the electrons have a 
much smaller mass compared to protons. This means that for a given pressure gradient of the gas the electrons tend 
to be accelerated much more than the ions. This leads in general to an electric field, which couples back positive 
and negative charges, of the form Et = — Vp e /ert e , where p e and n e are the electron pressure and number density, 
respectively. If such a thermally generated electric field has a curl, then by Faradays law of induction a magnetic field 
can grow. Taking p e = n e kT with T the electron temperature we have V x E^ = (ck/e)(Vn e /n e ) x VT. So E^ has a 
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curl only if the density and temperature gradients, are not parallel to each other. The resulting battery effect, known 
as the Biermann battery, was first proposed as a mechanism for thermal generation of stellar magnetic fields [0. 

The Biermann battery can also lead to the thermal generation of seed fields in cosmic ionisation fronts B . These 
ionisation fronts are produced when the first ultra violet photon sources, like quasars, turn on to ionise the intergalactic 
medium (IGM). The temperature gradient in a cosmic ionisation front is normal to the front. However, a component 
to the density gradient can arise in a different direction, if the ionisation front is sweeping across arbitrarily laid down 
density fluctuations, associated with protogalaxies/clusters since these in general have no correlation to the source of 
the ionising photons. The resulting thermally generated electric field has a curl, and magnetic fields on galactic scales 
can grow. They turn out to have a strength B ~ 3 x 10~ 20 G. This field by itself is far short of the observed microgauss 
strength fields in galaxies, but it can provide a seed field, coherent on galactic scales, for a dynamo. Indeed the whole 
of the IGM is seeded with small magnetic fields. This seed field may infact have the right symmetry properties for the 
galactic dynamo modes H . The Biermann battery has also been invoked to generate seed magnetic fields in galactic 
or proto-galactic environments Q . 

Larger seed magnetic fields than above can arise if we combine some form of dynamo action with the battery effect. 
For example, if stellar dynamos work efficiently, and some stars blow out as supernovae, then they can seed the 
interstellar medium, with significant magnetic fields. Alternatively galactic turbulence can itself lead to a small-scale 
dynamo (see below) and provide a larger seed for the large-scale galactic dynamo. There have also been attempts to 
invoke exotic physics in the early universe to produce primordial magnetic fields . (Infact a primordial field in the 
IGM, which redshifts to a present day value of ~ 1CP 9 G, and is correlated on Mpc scales, can significantly perturb 
structure formation, and cause detectable anisotropies in the cosmic microwave background radiation 0.) It is fair 
to say at present that most seed field generation mechanisms fall far short of producing large-scale correlated fields 
at the micro-gauss level. One does need some form of large-scale dynamo action. 



III. THE LARGE-SCALE GALACTIC DYNAMO 



Disk galaxies are differentially rotating systems. Also the magnetic flux is to a large extent frozen into the fluid. So 
any radial component of the magnetic field will be efficiently wound up and amplified to produce a toroidal component. 
But this results in only a linear amplification of the field. To obtain the observed galactic fields starting from small 
seed fields one should find a way to generate the radial component from the toroidal one. If this can be done, the 
field can grow exponentially and one has a dynamo. 

A mechanism to produce the radial field from the toroidal field was invented by Parker ||. The essential feature 
is to invoke the effects of cyclonic turbulence in the galactic gas. The interstellar medium (ISM) is assumed to be 
turbulent, due to for example the effect of supernovae randomly going off in different regions. In a rotating, stratified 
(in density and pressure) medium like a disk galaxy, such turbulence becomes cyclonic and acquires a net helicity. 
Helical motions of the galactic gas perpendicular to the disk can draw out the toroidal field into a loop which looks 
like a twisted f2. Such a loop is connected to a current and because of the twist this current has a component parallel 
to the original field. If the motions have a non-zero net helicity, then the random current components parallel to 
the field, add up coherently. A toroidal current can then result from the toroidal field. Hence, poloidal fields can be 
generated from toroidal ones. 

In quantitative terms, suppose the velocity field is the sum of a mean, large-scale velocity vo and a turbulent, 
stochastic velocity vy. The induction equation becomes a stochastic partial differential equation. The equation for 
various moments of B, can be derived in two ideal limits. First when R m << 1, and the distortions to the mean 
magnetic field are small, and second when R m >> 1, but the turbulence is assumed to have a delta function (or very 
small) correlation in time. For galaxies the latter idealisation may be more relevant. 

Let us split the magnetic field B = Bo + <5B, into a mean field Bo and a fluctuating component i5B. Here the mean 
is defined either as a spatial average over scales larger than the turbulent eddy scales or more correctly as an ensemble 
average. Assume the turbulence to be isotropic, homogeneous and helical. The action of the turbulent velocity field 
v-r, on the magnetic field, the (v x B) term, then leads to an extra contribution to the mean electric field of the form 
— cE = QfB — %V x B . Here a = —(1/3) J < vy(<).(V x vy(s)) > ds, depends on the helical part of the turbulent 
velocity correlation function, and r\ t — (1/3) J < vt(4).vt(s) > ds, called the turbulent diffusion depends on the 
non-helical part of the turbulence. Here the angular brackets <> denote an ensemble average, over the stochastic 
velocity field. The induction equation for the mean field, with the extra turbulent component of the mean electric 
field, is then given by 

<9B 

— - = V x (v x B + aB - fa + n T )V x B ) . (2) 
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This kinematic mean- field dynamo equation, can have exponentially growing solutions, which have been studied 
extensively in the literature fl. While the a-effect is crucial for regeneration of poloidal from toroidal fields, the 
turbulent diffusion turns out to be also essential for allowing changes in the mean field flux. Also in galaxies, 
differential rotation (the f2 effect) is dominant in producing toroidal from radial fields. The growth rates of the 
galactic 'a-f2 dynamo', are typically a few times the rotation time scales, of order 10 9 yr. Modulations of a, and t]t, 
due to enhanced turbulence along spiral arms, can also lead to bi-symmetric large-scale fields, correlated with the 
optical spirals |j) . 

Note that the mean field has a scale limited only by the size of the system, which can be much bigger than the 
scales associated with the turbulence. In this sense one has created order from chaos. One may be tempted to refer 
to this as an inverse cascade, a term which would suggest transfer of power from smaller to larger and larger scales. 
But reality is more subtle. All scales larger than the turbulent eddy scale can grow simultaneously due to the a-effect; 
but with larger scales growing slower. So the effect can be thought better as a long range interaction between the 
turbulent scales, and all larger scales. 

A physics comment is in order at this stage. When one considers the effect of turbulent fluid motions on a scalar 
field, like say smoke, one only gets a mean diffusion of smoke, associated with the random walking nature of turbulent 
fluid motions. But for 'frozen' magnetic fields the induction equation has terms which not only imply a body transport 
due to the random motions (v.VB), but also a term, B.Vv, which describes the generation of magnetic fields due 
to velocity shear. It is this qualitative difference between magnetic fields and smoke that leads to an alpha effect 
(provided also that the motions have a non-zero average helicity), over and above turbulent diffusion. As we see below, 
it also leads to the small-scale dynamo. Note that both a and r]T, depend crucially on the diffusive (random-walk) 
property of fluid motion. So if due to some reason (see below) the fluid motion becomes wavelike, then the above 
integrals average to zero, and the alpha effect and turbulent diffusion will be suppressed. 

In deriving the mean-field equation, the turbulent velocities have been assumed to be given, and unaffected by the 
Lorentz forces due to the magnetic field; at least not until the mean large-scale field builds up sufficiently. However 
this does not turn out to be valid due to the more rapid build up of magnetic noise compared to the mean field, a 
problem to which we now turn. 

IV. THE SMALL-SCALE DYNAMO AND MAGNETIC NOISE 

The problem with the kinematic mean field dynamo is that it is a myth. This is because, the same turbulence 
which contributes to the a-effect, the turbulent diffusion, and associated growth of mean fields, also leads to a more 
rapid growth of small-scale fields. 

In incompressible turbulence, fluid particles random- walk away from each other. This leads to stretching of the field 
lines attached to these particles, and an exponential increase of field strength. The stretching will also be accompanied 
by the field being squeezed into smaller and smaller volumes. Suppose one considers a flux tube, of length Z t , cross 
section A, density p and magnetic field B. Then flux freezing implies BA = constant, mass conservation pAl t = 
constant, and incompressibility p ~ constant. So as It increases due to random stretching, the the magnetic field 
B oc l t increases and the cross-section A cx decreases. 

If for the moment one ignores Lorentz forces, then the sqeezing into small volumes, stops only when diffusive scales 
are reached. Typically the field can be thought of as being in flux ropes, curved on the eddy scale, say L, and a 
thickness of order the diffusive scale, say (assuming only a single scale eddy is present). At this stage the energy 
input into the magnetic field due to random stretching would be comparable to the energy loss in diffusion. This 
gives vB/L ~ r\Bjr\, implying ~ L/Rm 2 - What happens further (whether growth or diffusion wins out), can only 
be decided by a more quantitative treatment of the problem. 

A rigorous analysis of small-scale field dynamics, was first worked out by Kazantsev ]l0| , and elaborated extensively 
by many authors |0], for the simple case, where the turbulence was assumed to have a delta function correlation 
in time. We shall mainly draw upon our own work in Ref. jl2| . It turns out that, the small-scale dynamo (SSD) 
can operate under fairly weak conditions; that the MRN associated with the turbulent motions be greater than a 
critical value R c ~ 100. In particular the fluctuating field tangled on a scale I, can grow on the turn over time scale 
of a turbulent eddy of scale I, with a growth rate T/ ~ vi/l, provided the MRN on that scale R m (l) — Vil/t) > R c . 
Here vi is the velocity associated with eddies of scale I. For Kolmogorov turbulence, since vi oc Z 1 / 3 , the growth rate 
T ; oc r 2 / 3 , and so increases with decreasing I. For galactic gas, with a significant neutral component, typically, even 
the eddies at the cut-off scale of the turbulence, say l c , have R m (l c ) >> R c . 

The spatial structure of the small-scale dynamo generated field is also of great interest. For this it is useful look at 
the behaviour of the magnetic correlation function, say w(r, t) =< <5B(x, <).<5B(y, t) >, where r = |x — y|. Here the 
averaging indicated by <>, is a double ensemble average over both the stochastic velocity and fluctuating magnetic 
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fields. From V.<5B = 0, one can show that the curve r 2 w{r) should enclose zero area. Since w(0) is neccessarily 
positive, as one goes to larger r, there must be some r ~ d say, when w(r) becomes negative. For such r, the 
fluctuating field at the origin, and at a separation d are pointing in 'opposite' directions on average. This can be 
intcpreted as saying that the field lines on average are curved on scale d. 

For Kolmogorov type turbulence, and if R m (l c ) >> R c , the fastest growing mode, has w(r,t) strongly peaked 

1 /2 

within the diffusive scale of the cut-off scale eddies, r = rd{l c ) = l c /Rm (l c ), changing sign at r ~ l c , and rapidly 
decaying for larger r/l c . One can interpret such a correlation function as implying that field is concentrated into 
ropes of thickness rj,(l c ) and curved on scales of order l c . For slower growing modes, with growth rate T/ ~ Vi/l, w(r) 
extends upto r ~ I, after which it decays exponentially. For these modes, the small-scale field can be thought of as 
being concentrated in rope-like structures with thickness of order the diffusive scale rd{l c ) and curved on a scale upto 
~ I. In general, at the kinematic stage, the growth rate of irreducible higher order correlations, increase with order, 
indicating that the field becomes highly intermittent in space. 

Note that the small-scale dynamo due to even the eddy at the outer scale of galactic turbulence, will lead to the 
exponetial growth of small-scale fields on a time r = L/v ~ 10 7 yr. (Here we have taken L ~ 100 pc, and a velocity 
scale v ~ 10 km s _1 .) This time is much shorter than the time scale ~ 10 9 yr for mean field growth. The magnetic 
field is then rapidly dominated by the fluctuating component, before the mean field has grown appreciably. If the 
energy in the small-scale component grows to equipartition with the turbulent energy density, the turbulence could 
become more wavelike 'Alfven' turbulence, than an eddy like fluid turbulence. So diffusive effects like the a and t]t, 
would get suppressed, and mean field growth stopped. How does the galaxy escape this predicament? 



V. SATURATION OF THE SMALL-SCALE DYNAMO 



To answer this question, it is crucial to find out how the small-scale dynamo saturates. We have concentrated on the 
possibility that the small scale field continues to be intermittent in space, when it saturates, as it was in the kinematic 
regime. That it saturates as a 'can of worms'; with peak fields being limited by non-linear effects to values of order 
or slightly larger than equipartition fields, but with most of the space having much smaller fields. Then the average 
energy density of the saturated small-scale dynamo generated field, may still be sub-equipartition, since it does not 
fill the volume. And the turbulence will remain eddy like, and preserve diffusive effects like a and rj t . We have given 
one explicit realisation of the above idea, in Ref. |12|, in the case of a galaxy dominated by neutral particles. 

In partially ionised plasma, the Lorentz force acts on ions, which are only coupled to neutrals through collisions. 
This leads to a 'ambipolar' drift of ions (and hence the field) with respect to neutrals. With v in Eq.(|l|) replaced by 
the neutral fluid velocity, the effective diffusivity changes to rf e ff = < <5B 2 > /(67rp,fj n ). Here pi is the ion density 
and Vi n is the neutral-ion collision frequency. So, as the energy density in the fluctuating field increases, the effective 
magnetic Reynolds number, for fluid motion on any scale of the turbulence say R a mu{l) — Vll/Veff, decreases. If 
RamU could decrease to R c , this itself will lead to a SSD saturation, but for conditions appropriate to galactic gas, 
RamU remains much larger than R c . So as the the small-scale field grows in strength, it continues to be concentrated 
into thin ropy structures, as in the kinematic regime. These flux ropes are curved on the turbulent eddy scales, while 
their thickness is now set by the diffusive scale determined by the effective ambipolar diffusion. 

Other restraining effects have to then limit the SSD. The first of these is due to the growing magnetic tension 
associated with the curved flux ropes. It acts to straighten out a flux rope, at a rate determined by equating the 
tension force and frictional drag. Frictional drag also damps the magnetic energy associated with the wrinkle in the 
rope. Further, small-scale flux loops can collapse and disappear, causing an irreversible sink of magnetic energy into 
heat. These non-local effects operate on the eddy turnover time scale, when the peak field in a flux rope, say B p , has 
grown to a few times the equipartion value. Their net effect is to make the random stretching needed for the SSD 
inefficient and hence saturate the SSD. As the field is in flux ropes which do not fill the volume, the average energy 
density in the saturated small-scale field is still sub-equipartition, and a and rjx are preserved. 

Note that B p has to grow to a larger and larger value, thinner the flux rope, for inefficient random stretching to 
operate. This is because tension is a volume force (oc r^lt) while drag acts on the surface of the rope and is oc r^lt- But 
B p cannot grow larger than {SirPext) 1 ^ 2 , where P ex t is the total pressure in the ISM. At the same time the thickness 
of flux ropes is larger, greater the ambipolar diffusion, or smaller the ion-density. So the SSD saturates as above, for a 
given P ex t, provided is less than a crictical value n\. In the ISM, if P ext is a factor F greater than the gas pressure 
(the gas assumed to have T ~ 10 4 K and a density n„), one gets n\ ~ lQ~ 2 cm~ 3 (v/l0krns~ 1 )~ 3 (n n /cm~ 3 ) 2 / 3 (F/2) 7 / 3 . 
So for a range of 'galactic' like parameters, this picture of small-scale dynamo saturation works. 

For larger ion densities n$ > n\ , the way the SSD saturates is not very clear. The peak field is still limited by the 
external pressure. But how the flux ropes behave in the post-kinematic stage is yet to be rigorously worked out. It is 
possible that, when one starts with weak fields, the field is first squeezed into small volumes until limited by magnetic 
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pressure. Subsequently, constructive folding of the field, may lead to fusing and thickening of the flux rope, while 
destructive folding may lead to reconnection, and dissipation. A phenomenological model fl3|| which incorporates 
this thickening of flux ropes as the field builds up, drives the SSD into saturation, when the rope thickness becomes 
of order LjR x J 2 , the peak field reaches equipartition levels, but with the average energy density of order R^ 1 of 
equipartition. 

Numerical simulations of dynamo action due to mirror-symmetric turbulence fli| or convection [ fl5[ have also 
hinted at a saturated state of SSD as described above; a magnetic field concentrated into flux ropes, occupying a 
small fraction of the fluid volume, having peak fields comparable or in excess of equipartition value but average 
magnetic energy density only about 10% of the kinetic energy density. Such simulations are however limited by the 
MRN they can achieve. There have also been MHD simulations of the SSD in fourier space, adopting some form 
of closure approximation, like the EDQNM approximation Jig] , or the DIA jf7|. They have also indicated that the 
small-scale field could saturate at sub-equipartition levels. However, why this happens and the relation of the fourier 
space, to the real space calculations which we have emphasised, is not at present clear. 

There ofcourse remains the possibility that the SSD in a fully ionised gas saturates only when the energy density 
in the noise grows comparable to that of the turbulence. However, the spatial structure of the small-scale field, is 
still likely to be highly intermittent. The large-scale dynamo action will then depend on how such a field responds to 
turbulent motions, especially whether the field can reconnect efficiently fLq| . If reconnection is efficient, then it may 
allow diffusive transport to still occur, through the forest of small-scale fields, rather like Tarzan, swinging from one 
rope to another crosses the jungle! Reconnection is an important issue, which deserves much more discussion than 
we have given (cf. [^9[). Another important issue which is just beginning to be addressed ]2C[], is the calculation of a 
(or t]t) in the presence of significant small-scale fields. 

In summary, the survival of the diffusive effects needed for large-scale dynamo action could depend crucially on 
whether the SSD generated fields can saturate at sub-equipartition levels. We feel that this may indeed be possible, if 
the noise saturates as a 'can of worms'. Note the SSD generated field is indeed spatially intermittent in the kinematic 
regime. So when one starts from weak fields, kinematic evolution operates for some time, and the small-scale field 
is driven to an intermittent state. The important question is to what extent it remains so in the non-linear regime, 
when it saturates. This SSD generated field will also provide a strong seed for the large-scale dynamo plj . Indeed, in 
a unified treatment of small and large-scale dynamos, large-scale correlations are produced, from small-scale fields in 
a way analogous to 'quantum mechanical tunnelling', of the stationary state of the SSD p3| . The helical, turbulent, 
a-fl dynamo still seems to be the best bet for explaining large-scale galactic fields. Ultimately the galactic dynamo 
is a non-linear dynamo; but a discussion of how it operates in the final saturated regime needs more thought. Our 
thoughts on galactic magnetism have yielded interesting results but much remains to be done. 
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